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In transportation networks, a spontaneous jamming transition is often observed, e.g. in urban 
road networks and airport networks. Because of this instability, flow distribution is significantly 
imbalanced on a macroscopic level. To mitigate the congestion, we consider a simple control method, 
in which congested nodes are closed temporarily, and investigate how it influences the overall system. 
Depending on the timing of the node closure and opening, and congestion level of a network, the 
system displays three different phases: free-flow phase, controlled phase, and deadlock phase. We 
show that when the system is in the controlled phase, the average flow is significantly improved, 
whereas when in the deadlock phase, the flow drops to zero. We study how the control method 
increases the network flow and obtain their transition boundary analytically. 


I. INTRODUCTION 

Transportation on networks plays an essentially pre¬ 
requisite role in our social activities. Complex network 
structures enable goods supply [J-Q, information com¬ 
munication 0,1, and human mobility mi in a vast 
human society. However, as a drawback of this efficiency, 
the structural complexity often causes breakdowns of the 
systems p~7l420| . In particular, traffic networks (e.g. air 
traffic networks [21] and road traffic networks Q) dis¬ 
play a macroscopic jamming transition when the system 
congestion exceeds a certain level. It is essentially dif¬ 
ferent from the classically known (microscopic) jamming 
phenomenon [22l . [23j and appears only when subsystems 
(e.g., roads or airports) are connected with each other. 
This macroscopic jamming is believed to explain (at least 
partially) severe jamming phenomena in reality [f||2l|. A 
simple guaranteed solution to the problem is to increase 
the effective capacity of the system, which is, however, 
not always feasible in our society confronted by rapidly 
increasing traffic demand. Hence, some external control 
method is necessary to mitigate the congestion. 

In such complex systems, understanding how auto¬ 
matic controls affect their macroscopic features still re¬ 
mains a grand challenge; nevertheless we have to deal 
with such exigent problems in reality. To facilitate the 
understanding of network traffic under control, innumer¬ 
able studies have been conducted (see e.g., Ref. |24j 
and references therein). However, because until very re¬ 
cently the issue of the macroscopic jamming transition 
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on networks has not been explicitly discussed, few stud¬ 
ies have focused on the outcome of control to mitigate 
such jamming. In this paper, we consider simple controls, 
which close nodes adaptively by temporarily removing 
links leading to congested nodes and reconnecting them 
when the congestion is resolved. This type of simple on- 
off controls is widely adopte d, e.g. i n traffic light controls 
[2~il IT! and ramp metering [24l. l27| for vehicular traffic, 
and the ground delay programs [28| for air traffic [29k 
[Ml . We show that this control method could cause a sud¬ 
den breakdown, which is characterized by a discontinuous 
phase transition to the deadlock state, whereas it is partly 
effective for alleviating congestion in a certain parameter 
regime. This deadlock transition stems from the control 
we adopt, and is a different phenomenon from previously 
known deadlock phenomena observed e.g. in the BML 
model [Hil. By disclosing the dynamics caused by 
the control, we discuss its effectiveness and criteria for 
designing a control law on transportation systems. 


II. MODEL 

We begin by reviewing the original model 0-H, [2l| of 
network transportation without controlling rules. It is 
known that flow of spatially exclusive objects often non- 
linearly depends on the density of the objects [H, HHHH ■ 
The flow-density relationship consists of two parts; i.e., 
free-flow and jamming regions, which are characterized 
by the increase and decrease of flow according to density, 
respectively. We consider density distribution on a net¬ 
work comprised of N nodes that are labeled by * (, j and 
l)= 1, • • • N. For each density value in node i, pi, flow on 
a link from the node is determined by the function value, 
J(pi). For simplicity, we assume a piecewise linear flow- 
density relationship [Fig. Oja)], J(pi) = min{pj, 1 — p^}. 



2 


(a) (b) 



(c) 



FIG. 1. (a) Flow-density relationship and control points, (b) Inflow to and outflow from a node in the network, (c) 

Disconnection of a link leading to a congested node. 




FIG. 2. (a) Queuing rule, (b) Detouring rule. 


Each node sends and receives density through the links 
[Fig. QJb)] defined by an adjacency matrix A, whose el¬ 
ement Aj,j is 1 if there is a link from node i to j, and 
0 otherwise. Hence, time development of the density in 
node i is written as 


dpi 

dt 


N 


N 


y ] Ajjj(pj) y ] Ajjj(pi). 

j =i t=i 


(i) 


In this paper, since we are interested in the nature 
of the system under control, we ignore individual prop¬ 
erties of nodes; the flow-density relationship, capacity 
and degree of nodes. Moreover, for simplicity directed 
random networks with homogeneous in- and out-degrees 
ki n = k ou t = k are used. 

It is known that this model [Eq. (©] exhibits spon¬ 
taneous destabilization for large densities [M3. When 
the density in a node is in the jamming regime (p > i), 
a small increase in density leads to a decrease of outflow, 
which causes further congestion. This effect destabilizes 
uniform flow in networks (for a more mathematical ex¬ 
planation, see Appendix 0. Hence, we consider a con¬ 
trol method as follows. Since when the density in the 
system is high, density concentrates in congested nodes, 
the method must disperse it in some way. 

As the simplest control, we here consider the following 
operations: (i) when the density in node i exceeds p c i, 
links leading to the node are disconnected [Fig. [2(c)], 
i.e., Aji = 1 —*■ Aji = 0 (j ^ i) (closed); (ii) when the 
density in a closed node drops below p op , the incoming 
links are all reconnected [see also Fig. [2a)]. The flow 
that is supposed to travel through the disconnected links 


must stay in the present node or detour using other links. 
In this paper, both rules are studied separately. The first 
rule is the queuing rule [Fig. [2a)], wherein the canceled 
flow stays in the departure node, resulting in no change 
in the system equation for open nodes Q, whereas for 
closed nodes the inflow term vanishes: 


dpi _ J AjiJ(pj) J2j=i A ijJ{Pi) (* ■ open), 

dt | - J2f=i AijJ(pi). (i : closed) 

( 2 ) 

The second rule is the detouring rule [Fig. (2b)], which 
redistributes the canceled flow equally to open links. 
Thus, flux to an open node is increased by this additional 
flow: 


dpi 

dt 


E N kj(pj) 

i =1 iAWO 
-kJ(pi) 


kJ(pi) 


(i : open), 

(i : closed), 

( 3 ) 


where o(l) represents the state of node l, i.e., open 
(o(l) = 1) and closed (o(l) = 0) states. For closed nodes, 
the equation is not influenced by the state of neighbor¬ 
ing nodes if at least one outgoing link is open. When the 
destination nodes are all closed, the outflow from a node 
vanishes, regardless of its state. 


A. Equivalence of the two control rules 

In a mean-field treatment (large N and k; see Ap¬ 
pendix the system equations of the two rules are 
equivalent through a transformation, — = t^, where 

x = is the ratio of open nodes in the thermo¬ 

dynamic limit. Hence, both rules result in the identical 
density distribution in the mean-field limit. 


B. Parameter settings 

In this paper we restrict ourselves to considering the 
case of pci = 0.75, and focus on the system’s depen- 

dence on p op and the average density, p = ■ 

We conducted simulations for N = 100 and k = 10 to 
compare the results with the predictions obtained by a 
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FIG. 3. Snapshots of (a) the free-flow phase (p = 0.4), (b) the controlled phase (p = 0.55) and (c) the deadlock phase 
(p = 0.65). Each disk represents a node, whose radius and color correspond to the density value. When a node is closed, the 
node and disconnected links are expressed by a magenta disk and gray lines, respectively. For display purpose, the networks 
are generated for IV = 50 and k = 3. 


mean-field theory. Note that the arguments in this pa¬ 
per are not affected significantly by these specifications. 
(See also Appendix [Cj) The system equations were in¬ 
tegrated by using the Runge-Kutta fourth-order method 
with dt = 0.0001. 


III. RESULTS 

A. Free-flow, controlled and deadlock phases 

First, we focus on the case of p op = 0.5. Depending 
on the value of the average density, the system exhibits 
three different phases, i.e., the free-flow, controlled, and 
deadlock phases (Fig. [3]). These phases are observed in 
the system both under queuing and detouring operations. 
When the average density is small, congested nodes do 
not appear and uniformly distributed flow is achieved 
[Fig. [3Ka)]. As the average density increases, congested 
nodes appear and our control rules apply. Because these 
operations increase the load in neighboring nodes, con¬ 
gested nodes appear on different locations while the oper¬ 
ated node is recovering from the congested state. Thus, a 
certain ratio of nodes are almost constantly closed; such 
a state is referred to as the controlled phase. When we 
further increase the average density, all the nodes be¬ 
come closed. In this state, nodes cannot accept nor send 
density, and thus flow in the system vanishes. Since the 
node capacity is temporarily limited during node closure, 
the operations effectively reduce the total capacity of the 
system. 

Next, we discuss the phase transitions for various p op 
values. Fig. [4] shows the phase diagram in the p-p 0 p 
plane. The simulation results were obtained by observing 
the state of the system after 10 6 simulation time steps. 
The initial density distribution in the system was set uni¬ 


formly, but perturbed with randomly selected 10 closed 
nodes (p = p c \). When p op is small, the controlled phase 
is not observed. In this region, closure time is too large, 
and thus it adversely affects the system by utilizing less 
of the capacity of the network. For large p op values, the 
controlled phase exists when the average density is mod¬ 
erately large. Here the control operations function with¬ 
out causing a breakdown. As a result, interestingly, the 
phase diagram has a triple point, which reflects the limit 
point of controllability. Because large p op values lead 
to small closure time of controlled nodes, the controlled 
phase region expands for such p op values. 


In a practical context, the transition line from the con¬ 
trolled phase to the deadlock phase is important to pre¬ 
vent a breakdown. Since the two rules are equivalent in 
the mean-field limit, we analyze the detouring rule, in 
which the time development of density in a closed node 
is known [Eq.©]. First, closed nodes in the controlled 
phase are focused on. We assume that the probability 
distribution of density in the time development, Q{p)dp , 
is proportional to the time in which the density is in the 
interval [p, p + dp] in the closed state, i.e., Q(p) oc 


dt 

dp 


From Eq. ©, we obtain 


dt 

dp 


= k j( p ) ■ Thus, the distri¬ 


bution is finally written as 


Q(p) = T c7io P 


l 

Wp)' 


( 4 ) 


where the normalization constant, T c i_y op , corresponds to 
the time from closure to the next opening of a controlled 
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B. Macroscopic flow-density relationship 
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if log 


1 

4(1 Pci) Pop 
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(Pop < 2 )> 
(Pop ^2)’ 

(Pop < 2 )j 
(Pop A 2 )- 


(5) 


Using this distribution, the average density in closed 
nodes is calculated as 


Pclosed — 


pQ{p)dp 


' Pop 


Pcl+Pop+log2(l-p„l)-l 
log4(l-pcl)p op 
1 Pci - Pop 

logi^SE’ 

6 i-Prl 


(Pop < 2 ) ; 
(Pop > 2 )‘ 


( 6 ) 


At the transition point when the deadlock occurs (all the 
nodes become closed), this value should coincide with the 
overall average density, i.e., p ~ p c iosed- This expression 
is shown in Fig. [H and gives an accurate prediction when 
Pop is large. Because this approximate analysis assumes 
that the number of closed nodes continuously increases 
from the controlled phase to the deadlock phase with in¬ 
creasing average density, the equation cannot describe 
the transition from the free-flow phase to the deadlock 
phase for small p op . The agreement between this equa¬ 
tion and the simulation results implies that the closure 
operations can be regarded to occur almost at random 
without causing synchronization. 



FIG. 4. Phase diagram of the system. The closed and open 
circles represent transition points obtained by simulations. 
The black bold line is the theoretical prediction, Eq. ©. The 
simulation results were obtained by examining the presence of 
a closed node after a sufficiently long time (f = 100). At t = 0, 
10 nodes were selected randomly to be closed with p — 0.75. 
With this initial condition, the system does not accidentally 
fall into the free flow (absorbing) state for system parameters 
in the other two phases. 


As the most important performance index of the sys¬ 
tem, we study a macroscopic relationship between the 
average flow per link and the average density, which is 

defined by J = ^'' t=1 • Figures 0D(a) and [51(b) 

show the macroscopic flow-density relationship for the 
queuing rule and detouring rule, respectively. It is conve¬ 
nient for understanding the results to consider the base 
line case p op = p c i = 0.75, wherein nodes reject inflow 
only when their (prescribed permissible) capacity (p c i) is 
fully occupied. Because of the system’s intrinsic insta¬ 
bility, the uniform density distribution is not achieved 
in the density region p > in which the macroscopic 
flow decreases. If the macroscopic flow is larger than 
that in the base line case, the controls are effective. In 
some cases, we can partly see the effectiveness in the con¬ 
trolled phase; however, when the average density further 
increases, the macroscopic flow sharply drops to zero, 
which corresponds to the phase transition to the dead¬ 
lock phase. When p op is small, the closure time is large, 
and thus the system falls into the deadlock phase even 
with a low average density, as can be seen in the phase 
diagram (Fig. 0j). In both rules, the average flow en¬ 
tirely increases with increasing p op when p op < i. As we 
further increase p op , the decrease in J in the controlled 
phase becomes gentle, and it can be achieved for a larger 
average density. Therefore, from a practical perspective, 
the most desirable p op depends on p (and p c \). 


In contrast with the queuing rule, the detouring rule 
does not stop flow even if the next node is closed by re¬ 
distributing it to other open nodes. Hence the average 
flow is larger than that operated under the queuing rule. 
This relationship is also verified by a scaling transforma¬ 
tion between them; = J<j with x < 1 (Eq. IB3I) . 


In the detouring rule, high average flow is achieved in 
the controlled phase. This can be understood by consid¬ 
ering the density distributions of nodes, P(p). Figure [G] 
shows the density distributions for p op = p c \ = 0.75 (the 
base line case) and p op = 0.5. In the base line case, the 
distribution is bimodal, whose peaks are approximately 
at p = 0.75 and 0.4. In this state, once a node becomes 
congested (p ~ 0.75), it cannot escape from the state be¬ 
cause the neighboring nodes with small density (p ~ 0.4) 
send more flow to the node than the outflow from the 
congested node. Thus, this density separation occurs. 
In contrast, for p op = 0.5, the density distribution is 
unimodal and the density separation is effectively sup¬ 
pressed by the closure operation that discharges density 
in congested nodes. Since the flow between nodes, J(p), 
takes the highest value in the middle density, collecting 
the density into the middle contributes to increasing the 
average flow that is given by integrating P(p)J(p). 
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FIG. 5. Macroscopic flow-density relationship for (a) the 
queuing rule and (b) the detouring rule. Note that, here 
the initial density values are set as pi = p + Si, where Si 
is a small random value drawn from a uniform distribution 
[-0.005,0.005]. 
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does not accept inflow). For this reason the phase dia¬ 
gram has a triple point. 

As expected, the random network structure in our 
study enabled us to analyze and understand the system 
clearly. Although this paper studies a limited range of 
networks, we have confirmed that the reported character¬ 
istics are not significantly influenced by inhomogeneous 
degree distributions (in random networks and an empiri¬ 
cal airport network in Ref. 2lj). Different network struc¬ 
tures would lead to changes in phase transition points 
(i.e., the effective traffic capacity of the system under 
control), which should be investigated in detail in future 
works. We believe that they will trigger rich varieties 
of studies such as tolerance of transportation networks 
[35, 136}], selective protection of nodes [3^-l40j. and the 
effects of spacial structure in a link |41h43I ]. 
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Appendix A: Stability of stationary flow 


FIG. 6. Density distribution (solid lines) and flow distribution 
(red broken lines) in the controlled phase (p = 0.6). 


IV. DISCUSSIONS 

A simple control method for mitigating congestion on 
network transport was investigated. We found complex 
behavior such as the deadlock phase transition, which 
is crucially important for implementing these rules to 
reality. The control impedes density separation among 
nodes, thus enhancing macroscopic performance of the 
system. However, to prevent deadlock situations, one 
can see the trade-off relationship between higher flow and 
small deadlock density. These intuitively understandable 
findings provide a useful overview for designing such con¬ 
trol systems. 

The phase diagram obtained in this paper is interest¬ 
ing also from a physical perspective. The deadlock phe¬ 
nomenon is similar to gridlocks observed in urban traffic 
m or a more simplified particle model, the BML model 
32H34j], which are caused by physical exclusion of vehi¬ 
cles (particles). In contrast, the deadlock in our model is 
caused by virtual exclusion of the controlling low (that 


Consider the state in which density is uniformly dis¬ 
tributed in the network, i.e., pt = p (i = 1, ■ • • , N). By 
linearizing the system equation with small perturbations 
(ei, ■ • ■ , ey) T = e added to this state, we obtain 

j t e = J'(p)Me, (Al) 

where matrix M is defined by Mjj = Aji [i ^ j) and 
Mu = —k. If all the eigenvalues of the matrix J'(p)M 
are not positive, the state is stable. Otherwise, a small 
perturbation grows to cause disruptions of flow. We show 
that the eigenvalues of M are not positive, and thus the 
stability of the system is determined only by J'(p). 

Let I be a matrix defined by L = M + kl, whose 
eigenvalues are Ai, • • • , Ay. Here I is the identity ma¬ 
trix. Since each element of this matrix is nonnegative, 
and the matrix is irreducible [Hj|, one can apply the 
Perron-Frobenius theorem f~15} to this matrix. Thus, we 
find that the spectrum radius of the matrix is k, i.e., 

| Ai|, • ■ ■ , | Ay | < k. Therefore the eigenvalues of the orig¬ 
inal matrix M, pi = Aj — k, are not positive. 

In conclusion, the system state is stable if J'(p) > 0 
and unstable if J'{p) < 0. 
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Appendix B: Relationship between the queuing and 
detouring rules 


Appendix C: Influence of degree k on system 
dynamics 


Here we apply a mean-field approximation to the links 
of each node as follows. We assume that k directed links 
starting from each node are composed of xk finks ending 
in open nodes and (1 — x)k links ending in closed nodes. 
Similarly, k directed links ending in each node are al¬ 
ways composed of xk links starting from open nodes and 
(1 — x)k links starting from closed nodes. In short, con¬ 
nections between links are assumed to be completely ran¬ 
domized, using x being the ratio of open nodes. Under 
this assumption, the time development of pi under the 
queuing rule obeys 


dpi. _ j -kxJiyPi) (i: closed), 

dtq {J2jLi Aj i J(p j ) - kxJ(pi) (i: open), 1 

and under the detouring rule, it obeys 

dfh ( -kJ(pi) (i: closed), 

dtd ~ \ T,f=i A ji J (Pj)/x ~ kJ{pi) (i: open), 


where subscripts q and d denote the queuing and detour¬ 
ing rules, respectively. It is straightforward from Eqs. 
(IB 1 1) and (IB2I) that in the mean-field approximation, the 
queuing and detouring rules cause the identical time de¬ 
velopment with a scaling transformation, xt q = td■ Cor¬ 
respondingly, the following relationship holds for average 
flow rates under the two rules. 


Figure H] displays the flow-density relationship for dif¬ 
ferent degree values, k. Difference in k does not influence 
qualitative characteristics of the system. As k increases 
the deadlock transition sharpens. When k is small, flow 
is often completely blocked because of the small number 
of routes from a node. As a result, the average flux is 
suppressed, but at the same time, due to the restriction 
of flow, the deadlock transition occurs at higher density. 



(B3) 


Figure [7] exhibits the comparison between and x q: 
which should be the same value if the two rules are equiv¬ 
alent. The agreement in the figure shows this equivalence 
is valid when the degree, fc, is large. 




FIG. 8. Flow-density relationship for different degree, k. Here 
we set Pop = 0.5. Other conditions are identical to those of 
Fig. 4(b). 


FIG. 7. Comparison between -j 1 and x q . Initial conditions 
are identical to those of Fig. 4. 
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